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We demonastrate the experimental technique for generating a spatial single-mode broadband 
biphoton field. The method is based on a dispersive optical element which precisely tailors the 
structure of the type-I SPDC frequency angular spectrum in order to shift different spectral com¬ 
ponents to a single angular mode. Spatial mode filtering is realized by coupling biphotons into a 
single-mode optical hber. 
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I. INTRODUCTION 

Progress in preparation and manipulation of quantum 
states of light is boosted by development of quantum 
communication and quantum computation. Spontaneous 
parametric down-conversion (SPDC) |T] is the most ef¬ 
ficient and widespread source of correlated photon pairs 
(biphotons). Phenomenologically, one may describe the 
SPDC process as a spontaneous decay of pump photon 
p traveling through a nonlinear medium into a pair 
of photons conventionally referred to as the signal s and 
idler i. Emerged photon pairs obey the energy conserva¬ 
tion law and phase matching relation: 

U!p=UJs+l^i, ( 1 ) 

kp = kg k-i ( 2 ) 

Here, ujp^g^i refer to the pump, signal and idler frequen¬ 
cies respectively, kp^g^i - corresponding wave vectors, and 
A is a phase mismatch responsible for the decay proba¬ 
bility [T]. In the present work we will concentrate on in¬ 
vestigating the frequency and angular characteristics of 
biphotons generated in the SPDC process. A biphoton 
quantum state may be expressed as 

1^) = F {ujg,eg,Wi,Bi)al {iOg,eg)a\ {uJi,9i) \vac), (3) 

where aj (wg, 9g) and a\ (wi, 9i) represent photon creation 
operators in frequency modes Ug^i and angular modes 9g^i, 
and F (wg, 9s,uji,9i) is a spectral amplitude of the bipho¬ 
ton. Relation between angular and frequency modes of a 
pair in monochromatic plane-wave pump approximation 
is given by: 

n{uji)ujisin9i = n{ujs)ujgsin9g, (4) 

where n(a;) is refractive index. Hence, the two-photon 
spectral amplitude F {ojg^9g,u}i,9i) reduces to a func¬ 
tion of two variables F{ujs,9g). The typical X-shape 
of the frequency-angular spectrum of the biphoton field 
\F (wg, 0s) P generated under collinear degenerate type-I 
phase-matching conditions is depicted in Figj^a. For ap¬ 
plications purposes a spatial single-mode biphoton field 
generation regime is preferential. A pump beam waist, 
which is typically IF > 50 /rm, defines photon pair emis¬ 
sion cone A9 — . This relation implies a narrow angu¬ 

lar emission range in the single-mode generation regime. 

Spectral broadening of the single spatial mode bipho¬ 
ton field remains a really significant puzzle among a 


wide list of biphoton spectrum control problems. The 
broadband biphotons exhibit ultra-short temporal corre¬ 
lations m- Their correlation time At ^ 1/A^. Ultra- 
correlated biphotons can be utilized in metrological appli¬ 
cations such as distant clock synchronization [3], quan¬ 
tum optical coherence tomography (QOCT) O |6] and 
quantum interferometric optical lithography BIS]. 

Besides, spectral broadening increases the degree of 
biphoton entanglement [5]> which can be exploited in 
quantum communication tasks. A quantum information 
can be directly encoded in frequency bins nsmn or the 
broad spectrum can be used for a wavelength-multiplexed 
entanglement distribution HMSj. 

Also biphoton spectrum broadening increases a maxi¬ 
mum pair generation rate Rmax ~ Av HZ] enabling an 
extremely high-speed quantum communication. 

Typically a biphoton spectral bandwidth is limited by 
an interaction medium dispersion A(a;) < 27r/L, where 
L is an interaction length. Hence, additional biphoton 
spectral broadening methods should be developed. Ho¬ 
mogeneous spectral broadening can be achieved by using 
the special nonlinear media with zero dispersion HSHU] 
{At = 5 fs [IH] and Av = 100 THz jZO] were obtained 
for type-I SPDC). Reducing nonlinear interaction length, 
one can generate 100 THz bandwidth type-H SPDC [22] 
and 150 THz type-I (23], but this method essentially de¬ 
creases an intensity of biphoton radiation. Pump beam 
focusing |24H26| and pump spectrum broadening [271125] 
also enable the broadband biphoton generation with a 
bandwidth up to 83 THz [25] . One more method of ho¬ 
mogeneous SPDC broadening is based on the modifying 
a crystal’s effective dispersion by means of a pump pulse- 
front tilt combined with a spatial walk-off [22]. A 44 THz 
bandwidth type-H SPDC spectrum was measured and a 
197 THz bandwidth (At = 6.4 fs) spectrum was theoret¬ 
ically predicted for type-I SPDC. 

Biphoton generation in inhomogeneous nonlinear me¬ 
dia (a set of crystals EHI, chirped crystals BE], crystals 
with a spatially modulated refractive index [32ff51] ) leads 
to inhomogeneous spectral broadening. Up to 154 THz 
bandwidth of the SPDC spectrum was obtained |34| . 
but such inhomogeneous broadening inevitably leads to 
a Fourier-unlimited field. 

All listed methods aim to increase parameter range (an 
area on the uj-9 coordinate plane) satisfying the phase¬ 
matching conditions which inevitably leads to a propor¬ 
tional drop in the spectral intensity. The typical view of 
the frequency-angular spectrum for broadened (namely 
by reducing a crystal length m) type-I SPDC is shown 
in Fig{T]b. 

It was demonstrated in BE] that the up-conversion 
process evinces an ultra-short correlation time (Ar = 
4.4 fs) provided that biphotons were emitted into a wide 
angular range and then refocused onto the second non¬ 
linear medium. However, collecting a wide angular range 
A9 results in essentially a multimode biphoton genera¬ 
tion regime, which is highly undesirable for long-distance 
(fiber or free-space) optical communication, microscopy 
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Figure 1. (Color online) The type-I collinear degenerate frequency-angular SPDC spectra, calculated for the 2 mm BBO crystal 
(zero on the frequency scale corresponds to the central frequency Vo = Wp/dvr). (a) The initial spectrum, (b) The typical view 
of the spectrum for usual broadening techniques, (c) The spectrum, transformed by an angular dispersion. The frequency 
bandwidth Ai/ corresponds to the collection angular range A0. 


and interferometric applications. Hence, the spatial 
single-mode broadband biphoton generation technique is 
a subject of endeavor for high-dimensional quantum-state 
engineering problems. 

We have to note that the spatial single-mode biphoton 
generation problem (leaving spectral properties aside) 
has been thoroughly studied (for example |36H39j ). Gen¬ 
erally, the solution is reduced to optimization of biphoton 
coupling conditions to a single-mode fiber. This implies 
maximization of single photon and coincidence counting 
rates and their ratio (the single to coincidence counting 
rate ratio is a key parameter in terms of pure single¬ 
photon state preparation). Our efforts were solely con¬ 
centrated on achieving the maximal spectral width of a 
spatial single-mode biphoton source. 

II. METHOD DESCRIPTION 

In spite of the fact that the total biphoton spectrum 
is broad both in frequency and angle, there is one-to- 
one correspondence between angles and frequencies due 
to phasematching. We propose to reconfigure the SPDC 
spectrum by means of an angular dispersive optical el¬ 
ement, shifting different spectral components to a sin¬ 
gle angular mode (Fig. [^c). This method doesn’t de¬ 
crease the spectral intensity of SPDC radiation and sig¬ 
nificantly increases the integral intensity in the target 
angular mode. 

The principal scheme of the experimental setup is de¬ 
picted in Fig. Different biphoton spectral compo¬ 
nents propagate in different angular modes. The lens 
system focuses emission onto an angular dispersive ele¬ 
ment (diffraction grating). Such a system transforms the 
SPDC angular spectrum to provide congruence with the 



Figure 2. (Color online) The principal scheme of the experi¬ 
mental setup 


grating dispersion curve. After diffracting from the grat¬ 
ing, the majority of spectral components correspond to 
the same angular modes and the radiation is coupled by 
the objective to a single-mode fiber. 

We have to point out that an analogous method was 
previously studied for optical parametric amplifier band¬ 
width broadening (for example |40fl44| ). However, it has 
never been utilized in biphoton generation applications. 
Nevertheless, biphoton spectral broadening using an an¬ 
gular dispersion module was demonstrated in |29| , where 
the pump pulse-front tilt was implemented by means of 
a diffraction grating. 


III. MODEL 

The numerical simulation of the SPDC process was 
evaluated based on the theoretical approach elaborated 
in [3H]. We consider only the photon pairs generated 
in given gaussian signal and idler modes, defined by a 
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Figure 3. (Color online) The pump, signal and idler gaussian 
modes. 


registration scheme: 


Ej {r) ~ gj (r) 


exp (ikjZj) exp 



+ y^) \ 


(5) 


where j = i, s,p. The pump field transverse profile is 
also assumed to be gaussian (Fig. 1^. On the one hand 
the biphoton emission rate 

Rt (x\F{uj,s,ujci,6LJ,e,,,e,,,W,,W,,Wp)\\ ( 6 ) 


depends on the mode overlap integral, and on the phase¬ 
matching conditions on the other. Here 0cs,^ci are the 
central angles of signal and idler beams, Wg, Wi and Wp 
are corresponding beam waists. The parameter 5uj cor¬ 
responds to spectral bandwidth of detection scheme. 

The diffraction grating dispersion relation yields the 
expression for binding angles 0cs,i and frequencies uJcs,i'- 

. 27rcL> 

sintlo - smytlcs.i = -, \‘) 

^cs,i 

where Oq is the incident angle of a gaussian mode con¬ 
jugated with the single-mode fiber fundamental mode 
(Fig.[i, D is the diffraction grating groove density and 
7 is the lens system magnification. The signal and idler 
beam waists Wg^i are related with the wavelengths As,i 
through 


fF,., = rWf (A,,,), (8) 

where F is the lens and condensing objective optical sys¬ 
tem magnification and Wj is the fiber fundamental mode 
beam waist, which is wavelength dependent: 

Wf = (9) 

^ ttN.A. ^ ^ 

Here N.A. is the value of the fiber numerical aperture. 

Substituting expressions 0.0-® to (|6|) one can ob¬ 
tain the frequency dependence of biphoton emission rate 
Rt (oJcs) governed by experimentally tuned parameters: 
7 , F and Wp. 



Figure 4. (Color online) Biphoton frequency spectra. The 
value u = 0 corresponds to the central frequency uip/in. Solid 
lines correspond to ideal calculation, dashed line corresponds 
to exact calculation,taking into account the diffraction effi¬ 
ciency, circles correspond to experimental data. Black thick 
lines and circles refer to transformed spectrum and red thin 
to initial one. 


IV. CALCULATION RESULTS 


Numerical calculations were performed for the nonlin¬ 
ear 2 mm long BBO crystal cut for a collinear degener¬ 
ate type-I phasematching, pumped by the 325 nm laser 
beam with a waist Wp = 34 p,m. Optimal SPDC cou¬ 
pling conditions demand shrinking the pump beam waist 
to smaller values but, taking into account that for given 
nonlinear media spatial walk-off is equal to 50 gm, fur¬ 
ther pump beam waist reduction won’t enhance a SPDC 
to single-mode fiber coupling efficiency. The diffraction 
grating groove density was chosen as 600/mm leading to 
7 = 1.05 for optimal overlap between the phasematch¬ 
ing curve and the diffraction grating dispersion curve. 
According to (38], to achieve the maximum collecting ef¬ 
ficiency F has to be selected in order to satisfy the con¬ 
dition 


Ws = Wi = V2Wp = 48 


^m. 


( 10 ) 


which corresponds to the angular bandwidth A0 = 0.5° 
in Fig. Considering the spectral dependence of Ws,i, 
F was smected to satisfy (10) only for the central wave¬ 


length (650 nm). Initial and transformed calculated co¬ 
incidence spectra are plotted in Fig.(straight curve). 

We used a diffraction grating blazed at 750 nm in our 
experiment. The first order diffraction conversion coeffi¬ 
cient at 750 nm was equal to 70% and drastically dropped 
down to 25% at 500 nm. We introduced diffraction ef¬ 
ficiency to our numerical model; the result is plotted by 
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Figure 5. (Color online) Biphoton correlation functions 
( |12[ | calculated from the spectra plotted in Fig. Solid lines 
correspond to ideal calculation, dashed line corresponds to 
exact calculation,taking into account the diffraction efficiency, 
circles correspond to experimental data. Black thick lines and 
circles refer to transformed spectrum and red thin to initial 
one. 


Figure 6. (Color online) The frequency bandwidth (111 de¬ 
pendence on the beam waist. Straight line - the ideal case, 
dashed line - taking into account the spectral dependence 
of the grating efficiency. Vertical line - the waist value 
VF = 48 jj,m used in the calculations and in the experiment. 


Table I. The spectral bandwidth Av ( |11[ ) and the second order 
correlation time At for the plots in Fig. and Fig. 



Initial 


Transformed 



Theory 

Exp. 

Theory 

Ideal Exact 

Exp. 

Au (THz) 

36 

39 

144 112 

99 

At (fs) 

20.3 

18.2 

3.7 4.1 

4.8 


a dashed curve in Fig. The values of the spectral 
bandwidths, determined as an integral characteristic: 


Ai/ = 


R[v)du/R{Q)^ 


( 11 ) 


are presented in Table |T] One may note that spectral 
dependence of diffraction efficiency leads to narrowing 
of the bandwidth of transformed spectrum from 144 to 
112 THz. 

Next we evaluated the corresponding second order cor¬ 
relation functions 


G^^^(t) cx 


/ 


2 




( 12 ) 


which are presented in Fig. The corresponding corre¬ 
lation times are shown in Table HI The correlation times 
have been calculated similarly to Ideally, 3.7 fs 

correlation time can be achieved. It corresponds to 1.7 
optical cycles. Non-uniform grating spectral efficiency 
leads to longer 4.1 fs correlation time. 


We calculated the spectral bandwidth (111 with re¬ 
spect to SPDC beam waists Ws^i (650 nm) considering 
the relation for the pump beam waist (lOl. The results 
are plotted in Fig. The vertical line corresponds to 
Ws^i (650 nm) = 48 ^m, included in our numerical model. 
Evidently, this value appears to be non-optimal consider¬ 
ing maximum spectral broadening. Increasing the beam 
waist up to 500 /rm leads to broadening of the spectral 
bandwidth up to 213 THz in an ideal case and up to 
154 THz taking into account the diffraction efficiency. 
On the other hand, larger beam waists imply a signifi¬ 
cant drop in coincidence count rate. According to |38| 
the coincidence count rate scales as oc W~^. 


For beam waists smaller than 100 pm our numerical 
model predicts a decrease in the spectral bandwidth of 
the coincidence count rate. This effect is associated with 
shrinking the interaction volume for a noncollinear gen¬ 
eration regime in a nonlinear medium between the pump, 
signal and idler gaussian beams, which leads to depletion 
of the spectral intensity of distant spectral components. 
At the same time, the collection angular range increases 
allowing the vertical branch of the angular-frequency 
spectrum (Fig.j^c) to enhance the central spectral com¬ 
ponents and, hence, reduce the spectral width. Here we 
have to note that our model doesn’t apply to a tight-focus 
scenario, when the waist length is shorter than the crys¬ 
tal length. Moreover, it doesn’t account for broadening 
of the SPDC spectrum due to pump focusing p41 - E5] . 

Spectrum bandwidth reduction due to the beam waist 
decreasing is associated with angular range Ad narrow¬ 
ing, where the phasematching curve overlaps the grating 
dispersion curve. Hence, the more precise overlap can be 
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achieved for the narrower spectral range. 


V. EXPERIMENT 

Fig. [^illustrates the scheme of the experimental setup. 
Laser radiation at 325 nm generated by the HeCd laser 
was focused by a lens LO (/o = 145 mm) on the 2 mm long 
BBO crystal cut for collinear degenerate type I phase¬ 
matching. A dispersive prism and filter FI blocked par¬ 
asite arc radiation. The pump beam was blocked with a 
UV Filter F2. Biphotons were focused through a double¬ 
lens system LI and L2 (/i = /2 = 40 mm) on the diffrac¬ 
tion grating, described in the previous section. The lens 
positions (noted in Fig.j^in mm) were selected to set the 
magnification of the system to 7 = 1.05. The numerical 
aperture of the lens system exceeded 0.17 which corre¬ 
sponds to ±9.5° angular range. It corresponds to the 
spectral mismatch interval from -150 to 460 THz for the 
given the phasematching curve (Fig. [^a). To measure 
an initial spectrum, the diffraction grating was replaced 
by a broadband mirror M. 

We used a 20X objective lens 01 to couple first-order 
diffracted light to a single-mode fiber SMF (Thorlabs 
SM600). This fiber retains single-mode propagation in 
the 600-800 nm band. For shorter wavelengths of 500- 
600 nm the fiber maintains a few-mode lossless propa¬ 
gation, which may affect the coupling efficiency in the 
specified spectral range. For our scheme the effect of a 
non-single-mode coupling is negligible. However, for the 
near-infrared range starting from 800 nm the propagation 
losses become significant. That is why we only used aim 
long fiber patch-cord to minimize propagation losses. 

Single-mode fiber output radiation was collimated by 
the objective 02 and split by a 50:50 nonpolarizing 
broadband beam splitter BS and coupled with objectives 
03 and 04 to multimode fibers MMFl and MMF2 con¬ 
nected with Single-photon counting modules D1 and D2 
based on Si avalanche photodiodes. One of the channels 
after the beam splitter included a spectrograph prism 
block, capable of i5iz = 2 -7 4 THz (wavelength depen¬ 
dent) spectral band filtering. Sufficiently wide spectral 
filtering enables both high coincidence count rate and 
sufficient spectral resolution. 

Experimental data took into account all spectral de¬ 
pendent parameters of the measurement part: detector 
quantum efficiency and wavelength dependent spectral 
filtering bandwidth. Thus, the key non-ideality of the 
preparation part (grating efficiency) was taken into ac¬ 
count in the above theoretical model. 


VI. EXPERIMENTAL RESULTS 

The measured transformed coincidence spectrum is 
plotted with thick (black) circles in Figj^ Thin (red) 



Preparation 


SMF 


Measurement 



Figure 7. (Color online) The scheme of the experimental 
setup. 


circles refer to the initial spectrum measured after plac¬ 
ing the broadband mirror instead of the diffraction grat¬ 
ing. All the plots are normalized against their maximal 
values. The measured spectral intensities at the central 
frequency for initial and transformed spectra were differ¬ 
ent, which can be explained by the distinction between 
the grating efficiency and mirror reflection coefficient. 


Agreement between the experimental data and numer¬ 
ical simulation (accounting for diffraction grating effi¬ 
ciency) is observed. The 99 THz spectral bandwidth 
(111 was obtained. The signal drop on the edges of 
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the spectral range is attributed to geometric lens aber¬ 
rations. Notably, the measured initial spectrum looks 
slightly broader than the calculated one. It can be ex¬ 
plained by the small broadening due to pump focusing 
], which wasn’t taken into account in the theory. 


The second-order correlation functions (121 calculated 
from the experimental data are also plotted with circles 
in Fig.[^ The correlation time related to the transformed 
spectrum is 4.8 fs. All spectral bandwidths Ai/ and cor¬ 
relation times 6t are presented in Tablein the "Exp." 
columns. 


The impact of the single-photon impurity is quantified 
by the coincidence to single count rates ratio: 1 — Rc/Rs. 
We experimentally estimated the ratio Rc/Rs (at the 
central frequency) to be 2.5%. The low biphoton intake 
is attributed to significant optical losses in the measure¬ 
ment part of the setup and an inefficient fiber coupling 
scheme. This is a trade-off between the maximal broad¬ 
ening and maximal efficiency. Thus, the maximum possi¬ 
ble value of Rc/Rs is 45 ± 5% assuming any other losses 
are negligible. The preparation part of the setup solely 
introduces optical losses reducing the Rs/Rc ratio down 
to 17 ± 2%. 

We measured coincidence count rates = 33 Hz and 
Rtr = 19 Hz at the central frequency (Sv = 3.2 THz) 
for initial and transformed biphoton fields to ensure 
the biphoton spectral power density invariance. The 
ideal scenario implies Rtr/Rin = 1- The measured 
ratio Rtr/Rin = 0.58 may be attributed to a signifi¬ 
cant difference between the diffraction grating efficiency 
and the broadband mirror reflectance at the given fre¬ 
quency. Moreover, a discrepancy between the measured 
value of Rtr/Rin ancl the corresponding broadening fac¬ 
tor Avin/Avtr = 0.39 endorses the fact that in our exper¬ 
iment the spectral intensity drop isn’t induced by bipho¬ 
ton spectral broadening. 


are required HZl iSJ iS]. The technique, demonstrated 
in our experiment, provides a Fourier-limited broadband 
biphoton source. Additional dispersive broadening of the 
correlation function may arise after passing through the 
lens system and single-mode fiber [3]. This effect can 
be avoided by substituting lenses with parabolic mirrors 
and realising the spatial mode filtering with a pinhole. 
Thus, our method (reversed in time) can also be applied 
for precise measurement of second-order correlation time 
by exploiting the sum-frequency generation scheme. The 
resolution of the such scheme is limited by the phase¬ 
matching bandwidth in the nonlinear crystal. 

However, our method exhibits significant optical losses, 
which leads to a considerable additional single-photon 
component in the state The fact that the fabrica¬ 
tion of efficient broadband diffraction gratings is barely 
feasible makes these kinds of losses inevitable. On the 
other hand, the losses emanate from non-optimal cou¬ 
pling biphotons to a single-mode fiber. The maximal 
coupling efficiency requires each biphoton spectral com¬ 
ponent to populate a single Gaussian mode. It can be 
realized by irradiating a nonlinear medium with a tightly 
focused pump beam. In our case the above condition is 
impossible to achieve for all frequency-angular compo¬ 
nents at the same time. 

It should be noted that the demonstrated method en¬ 
ables broadband biphoton generation in a single spatial 
mode, whereas applications based on correlation func¬ 
tion measurements demand photons of a pair to popu¬ 
late distinct spatial modes. This problem can be solved 
by a combination of two coherent single-mode biphoton 
sources m- 

Finally, the demonstrated technique significantly 
broadens the biphoton frequency spectrum without de¬ 
pleting the spectral intensity (and, moreover, with en¬ 
hancing the total intensity) by redistributing energy be¬ 
tween spatial biphoton modes. 


VII. DISCUSSION 

The demonstrated biphoton source engineering tech¬ 
nique possesses several distinctive advantages. First, the 
method provides an ultra-broadband biphoton genera¬ 
tion regime. The measured 100 THz bandwidth and 5 fs 
correlation time corresponds to values enabled by the 
homogeneous broadening techniques HH Ha I2E]- The 
optimization of our scheme enables bandwidth up to 
144 THz, which have been demonstrated by using the 
inhomogeneous |33] or extra-thin [23] nonlinear crystals. 

On the other hand, our method is limited by the 
spectral region, where phasematching and grating dis¬ 
persion curves coincide, whereas other spectral broaden¬ 
ing methods, based on utilizing inhomogeneous nonlin¬ 
ear media [231 [3nil31] . are completely free of this lim¬ 
itation. At the same time, the biphoton field, gener¬ 
ated in inhomogeneous nonlinear media, cannot be kept 
Fourier-limited, and additional compression techniques 


VIII. CONCLUSION 

We proposed the novel method of generating a spa¬ 
tial single-mode ultra-broadband biphoton field based 
on the frequency-angular spectrum transformation by 
means of the angular dispersion. The efficiency of this 
method was demonstrated both theoretically and exper¬ 
imentally and our calculations are in a significant agree¬ 
ment with our experimental results. The spectral band¬ 
width Av > 100 THz and the correlation time At < 5 fs 
can be reached. The key feature of our method - spectral 
intensity invariance — was also demonstrated. 
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